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ABSTRACT 

It is the purpose of this study to demonstrate the viability and economy of Design-of-Experi- 
ments methodologies to arrive at microscale secondary flow control array designs that maintain optimal inlet per- 
formance over a wide range of the mission variables and to explore how these statistical methods provide a better 
understanding of the management of flow in compact air vehicle inlets. These statistical design concepts were 
used to investigate the robustness properties of “low unit strength” micro-effector arrays. “Low unit strength” 
micro-effectors are micro- vanes set at very low angles-of-incidence with very long chord lengths. They were 
designed to influence the near wall inlet flow over an extended streamwise distance, and their advantage lies in 
low total pressure loss and high effectiveness in managing engine face distortion. The term robustness is used in 
this paper in the same sense as it is used in the industrial problem solving community It refers to minimizing the 
effects of the hard-to-control factors that influence the development of a product or process. In Robustness Engi- 
neering , the effects of the hard-to-control factors is often called “noise”, and the hard-to-control factors them- 
selves are referred to as the environmental variables or sometimes as the Taguchi “noise” variables. Hence Robust 
Optimization refers to minimizing the effects of the environmental or “noise” variables on the development 
(design) of a product or process. In the management of flow in compact inlets, the environmental or “noise” vari- 
ables can be identified with the mission variables. Therefore this paper formulates a statistical design methodol- 
ogy that minimizes the impact of variations in the mission variables on inlet performance and demonstrates that 
these statistical design concepts can lead to simpler inlet flo w management systems. 

INTRODUCTION 

The current development strategy for combat air -vehicles is directed towards reduction in the 
Life-Cycle Cost (LCC) with little or no compromise to air-vehicle performance and survivability. This strategy 
has been extended to the aircraft component level, in particular , the engine inlet diffuser system. One method to 
reduce inlet system LCC is to reduce its structural weight and volume. Consequently, advanced combat inlet con- 
figurations are being made more compact (or shorter) to achieve weight and volume (and LCC) reduction. How- 
ever, compact S-duct diffusers are characterized by high distortion and low pressure recovery, which are 
produced by extreme wall curvature and strong secondary flow gradients within the boundary layer. These char- 
acteristics are also further aggravated by maneuvering conditions. The requirement to highly integrate or embed 
the propulsion system into the airframe often leads to conformal inlet aperture shapes which do not lend them- 
selves to high aerodynamic performance. These configurations also present a challenging environment for both 
fan/compressor surge margin and aeromechanical vibration. Interest in High Cycle Fatigue (HCF) research by the 
US aerospace community has been spurred by discrepancies between the expected durability of engine compo- 
nents compared to that actually experienced in the field. Recognizing that inlet distortion is a forcing function for 
vibration in the fan components, methods for increasing HCF Life Expectancy can be combined with techniques 
for inlet recovery and engine face distortion management. Therefore, to enable acceptable performance levels in 
such advanced, compact inlet diffuser configurations, microscale secondary flow control (MSFC) methods are 

being developed to manage the recovery, distortion, and HCF aspects of distortion.®"® 
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One of the most difficult tasks in the design of a MSFC array for optimal inlet operation is arriv- 
ing at the geometric placement, arrangement, number, size and orientation of the effector devices within the inlet 
duct to achieve optimal performance.These effector devices can be either mechanical (positive displacement) or 
fluidic. This task is complicated not only by the large number of possible design variables available to the aerody- 
namicist but also by the number of decision parameters that are brought into the design process. By including the 
HCF effects in the inlet design process, the aerodynamicist has a total of seven individual response variables that 
measure various aspects of inlet performance. These include the inlet total pressure recovery, the inlet total pres- 
sure recovery distortion at the engine face and the first five Fourier harmonic 1/2-amplitudes contained in the 
engine face distortion pattern. Each of these responses must be maximized, minimized, constrained or uncon- 
strained while searching for the optimal combination of primary design variable values that satisfy the mission 
requirements. The design task is further complicated by the existence of hard-to-control factors that affect inlet 
performance, i.e. the mission variables. The mission variables that cause the off-design penalty are, for example, 
inlet throat Mach number (engine corrected weight flow), angle-of-incidence and angle-of-yaw. While the aero- 
dynamicist does not know how the pilot is ultimately going to fly the aircraft, it is known how the mission vari- 
ables affect inlet performance under wind tunnel conditions. Traditionally, tolerance or robustness to the mission 
variables was accomplished only after the parameter design was completed, usually by accepting whatever off- 
design performance was delivered by the newly designed inlet system. Numerical optimization procedures that 
have been successful with some aerodynamics problems give little assistance to designing robust inlets since they 
are point-design procedures, usually with only one decision parameter. However, there is a branch of statistical 
Design-of-Experiments (DOE) methodology which integrates both traditional Response Surface Methods (RSM) 
and Robust Optimization Concepts (ROC) into a single optimization procedure. It presents new potential for fur- 
ther reduction of total quality cost over the traditional design approach. 

Taguchi^ coined the term Robust Parameter Design to describe an approach to industrial prob- 
lem solving whereby the product variation is reduced by choosing levels of the control factors (design parame- 
ters) that make the product insensitive to the changes in the noise factors that represent sources of variations. 
These noise factors in industrial design are often the environmental variables such as temperature and humidity, 
properties of the material, and product aging. In some applications, the factors measure how the consumer uses or 
handles the product. In the aerodynamic design of inlet systems, there is an analogous situation to the industrial 
design problem. As mentioned above, the design of inlet systems is usually accomplished at the cruise condition 
(the on-design condition) while variations from the cruise condition are considered as an off-design penalty. 
Because the mission variables cause variation from on-design performance, they can be identified with the noise 
or environmental variables in the analogous industrial design problem. Likewise, how the pilot flies the aircraft 
can be identified with how the consumer uses or handles the product. In the industrial problem, researchers must 
be able to control the environmental variables in a laboratory environment, even though they cannot be controlled 
at the production level or in the field. Likewise, the aerodynamic researcher can indeed control the mission vari- 
ables in the wind tunnel environment, however these variables cannot be controlled in flight (in the field). By 
making the analogy between the industrial design problem and the aerodynamic design problem, Robust Optimi- 
zation Concepts developed for industrial problem solving can be adapted to the design of inlet systems^, and in 

particular, design of optimal robust microscale secondary flow control arrays®"® which can simplify the inlet 
flow management system. 


NOMENCLATURE 


AIP Aerodynamic Interface Plane 

c Effector Chord Length 

CCF Central Composite Face-Centered 

CFD Computational Fluid Dynamics 

D Engine Face Diameter 

DC60 Circumferential Distortion Descriptor 

DOE Design of Experiments 

h Effector Blade Height 

HCF High Cycle Fatigue 

Fk/2 k th Fourier Harmonic 1/2-Amplitude 

FM/2 Mean Fourier Harmonic 1/2-Amplitude 

L Inlet Diffuser Length 

LCC Life Cycle Costs 

MSFC microscale Secondary Flow Control 

Mt Inlet Throat Mach Number 

n Number of Effector Vanes per Band 

PFAVE Average Inlet Total Pressure at AIP 


R Inlet Radius 

R cl Centerline Radius 

R ef Engine Face Radius 

R thr Inlet Throat Radius 

ROC Robust Optimization Concepts 

Re Reynold Number per ft 

RSM Response Surface Methodology 

UAV Unmanned Air Vehicle 

UCAV Unmanned Combact Air Vehicle 

X cl Axial Distance 

Yj j Generalized Response Variable 

Ym,(x Response Variable Summed over Mt and a 

Z cl Centerline Offset Displacement 

a Inlet Angle-of-Incidence 

P Effector Vane Angle-of-incidence 

AZ cl Inlet Centerline Offset 

y Inlet Angle-of-Yaw 
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RESULTS AND DISCUSSION 


Baseline Flow in the Redesigned M2129 Inlet S-Duct 


The redesigned M2 129 inlet S-duct used in this study was considered similar to the original 
DERA/M2129 inlet S-duct defined by AGARD FDP Working Group 13 Test Case 3 using Lip No. 3 and For- 
ward Extension No. 2. This inlet design was first proposed by Willmer, Smith and Goldsmith, ^ and has been 
used extensively in the US and UK to explore inlet flow control array design. The centerline for the redesigned 
M2 129 inlet is given by the equation 


Z 


cl 



1 - cos 



( 1 ) 


the radius distribution measured normal to the inlet centerline is given by the expression 


Rrf ^thr^ _ ^ 
R ef ~ R thJ 


i 


41 1 — — 


+ 1 


( 2 ) 


where R thr = 2.5355 inches, R e j = 3.0 inches, L - 15.0 inches, and A Z cl = 5.7809 inches. The rede- 
sign of the M2 129 inlet was such that the new inlet matches the static pressure gradients normally found in typi- 
cal UAV or UCAV designs. Therefore, the new inlet is more compact than the original M2 129 inlet S-duct. As a 
consequence, supersonic flow will develop in this inlet when the inlet throat Mach number increases much above 
0.70. 


Inlet Flow Control Array Design 

To manage the flow in the redesigned M2 129 inlet S-duct, a single band array arrangement of 
microscale effectors was placed in the upstream section near the inlet throat. See Figures (1). These microscale 
effectors were micro-vanes, the largest height being about the average height of the momentum layer just down- 
stream of the inlet throat or about 2.0 mm. The purpose of these micro-vanes was to create a set of co-rotating 
vortices that would quickly merge to form a thin layer of secondary flow that will counter the formation of the 
passage vortex pair. Since the height of the vane effectors were limited to 2.0 mm, a single-band arrangement of 

micro- vanes set at 5.0° angle-of-incidence was chosen to investigate the enhancing effect of increasing the vane 
chord length on distortion management, i.e. allowing the array design to influence the inlet flow over an extended 
streamwise distance for a design advantage 

The DOE approach followed directly from the objectives of the study and are reflected in the lay- 
out of the design factors listed in Table (1). The factor variables were the number of vane effector units (n), the 
micro- vane effector height (h), the micro- vane chord length (c), the inlet throat Mach number (Mt), and the inlet 
angle-of-incidence (a). Strictly speaking, the inlet throat Mach number and angle-of-incidence are mission vari- 
ables and, therefore, the noise factors that belonged with the environmental variables, i.e. the outer array in a tra- 
ditional Taguchi-style Robust Parameter Design. However, in this study, the throat Mach number and inlet angle- 
of-incidence were combined into the statistical DOE matrix with the control factors. This is called a combined 

DOE matrix array, which allowed greater economy than the traditional Taguchi approach^. The robust nature of 
the throat Mach number and inlet angle-of-incidence was investigated during the analysis phase of the data. Table 
(2) shows the variables that were held constant during this study. They include the effector vane thickness (t), the 
geometric angle-of-incidence of the micro vanes ((3), the inlet operating total pressure (Pt) and temperature (Tt), 
and the inlet angle-of-yaw (y). Table (3) displays the response variables for this study. They include the inlet total 
pressure recovery (PFAVE), the engine face distortion (DC60), and the first five Fourier harmonic 1/2-amplitudes 
of engine face distortion (Fl/2, F2/2, F3/2, F4/2, and F5/2). 

The DOE strategy selected was a Central Composite Face-Centered (CCF) DOE. This strategy 
resulted in 27 unique CFD experimental cases that are shown in Table (4). This DOE construct is called a com- 
bined array format because it contains both the factor (design) variables and the environmental (mission) vari- 
ables. Notice that these DOE cases covered a substantial range of possible flow situations over a wide range of 

throat Mach numbers from 0.30 to 0.70, and angles-of-incidences from 0.0° to 20.0°. This particular DOE, like 
most DOE strategies, varied more than one factor at a time. Further, this layout of 27 cases permitted the estima- 
tion of both linear and curvilinear effects as well as two-factor interactive or synergistic effects among the DOE 
factors. This CCF DOE strategy is superior to the traditional approach where only changing one variable at a time 
does not permit the estimation of the two-factor interactions. It is also more economical at 27 runs than a full fac- 
torial approach where the number of experiments would be 3 5 or 243 separate CFD cases. It is also more econom- 
ical than a comparable Taguchi approach requiring 15 x 3 = 45 experiments. 

A graphical representation of the Central Composite Face-Centered DOE used in the study is 
presented in Figure (2). The DOE cases are represented in this figure by the circular symbols, where the symbol 
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locations on the cube signify its factor value. This DOE is called a composite DOE because the organization of 
cases is composed of a fractional factorial part and a quadratic part. The fractional factorial part of the DOE is 

composed of one-half of the 2 5 possible cases, i.e. 32 possible factorial cases, which are represented by the eight 
corner locations in each of the four corner-cubes in Figure (2). Because only half the number of possible factorial 
cases are actually used in this DOE format (circular symbols), the layout is called a 1/2-fractional of the full fac- 
torial and is composed of 2 5 " 1 cases, or 16 separate CFD runs. The remaining cases in Figure (2) are the quadratic 
part of the DOE. The quadratic cases allow for the evaluation of the curvilinear effects. All together, there are a 
total of 27 cases in a Central Composite Face-Centered DOE with five factor variables. Notice the balanced lay- 
out of cases in Figure (2). The factor variables are represented by the axes of the individual cubes, while the envi- 
ronmental (mission) variables are represented by the different cubes. This layout of cases represents the smallest 
number of CCF DOE cases that allows for the evaluation of linear and curvilinear effects as well as all two-factor 
interactive or synergistic effects. 

Each of the 27 cases in Table (4) was run with a Reynolds-averaged Navier-Stokes code^ that 
allowed for numerical simulation of micro-vane effectors without the need to physically embed the vane effectors 
within the CFD grid structure. However, for the present study the individual vanes were incorporated into the half 
cylindrical grid structure. These micro-vanes all had a thickness of 0.138 mm. See Table (2). The computational 
grid was developed such that it reasonably resolved the boundary layer development on both the suction and pres- 
sure surfaces of each micro-vane in the array. Because wall functions were used in the calculations, the grid reso- 
lution for the individual micro-vanes was simplified. However, the boundary layer along the micro-vane edges 
was assumed to be negligible, and therefore not resolved in the computational grid. The half cylindrical grid 
structure was composed of three blocks: an upstream block, an effector section containing the micro-vanes, and a 
downstream block. See Figures (1). The computational half-plane grid varied in total number of mesh points from 
about 950,000 to 1,150,000 depending on the micro- vane configuration. All CFD calculations were accomplished 
assuming half cylindrical symmetry. A two-equation k-1 turbulence model was used in this study. The model con- 
sists of transport equation for the turbulent kinetic energy and turbulent length scale. The model includes a near- 
wall model and compressible corrections for high speed flows. To introduce an angle-of-incidence (a-distur- 
bance) into the flow analysis, the condition was imposed that the initial station have an angle-of-incidence com- 
ponent that approximated the measured angle-of-incidence flow field^ 10 l 

Optimal Flow Control Over the Mission Variable Range 

To illustrate the potential of RSM and Robust Optimization Concepts to design and optimize 
MSFC arrays, three mission strategies were considered for the subject inlet, namely (1) Maximum Performance, 
(2) Maximum Engine Stability, and (3) Maximum HCF Life Expectancy. The Maximum Performance mission 
minimized the inlet total pressure losses, the Maximum Engine Stability mission minimized the engine face dis- 
tortion, while the Maximum HCF Life Expectancy mission minimized the mean of the first five Fourier harmonic 
amplitudes, i.e. “collectively” reduced all the harmonic 1/2-amplitudes of engine face distortion. Each of the mis- 
sion strategies was subject to a low engine face distortion constraint, i.e. DC60 < 0.10, which is a level acceptable 
for commercial engines, and a constraint on each individual Fourier harmonic 1/2-amplitudes: Fk/2 < 0.015, k = 
1,2... 5. For each of three mission strategies, i.e Maximum Performance, Maximum Engine Stability, and Maxi- 
mum HCF Life Expectancy mission, an “Optimal Robust” (open loop control) and an “Optimal Adaptive” 
(closed loop control) array were designed to operate over an inlet throat Mach number range from 0.30 to 0.70, 

and angle-of-incidence range from 0.0° to 20.0°. The “Optimal Robust” array arrived at a single MSFC array 
which operated optimally over the entire throat Mach number and angle-of-incident range (open loop control). 
The “Optimal Adaptive” array optimized all the design parameters at each throat Mach number and angle-of-inci- 
dence. Thus the “Optimal Adaptive” array would require a closed loop control system to sense a proper signal for 
each effector and modify that effector device, whether mechanical or fluidic, for optimal inlet performance. The 
inlet throat Mach number and angle-of incidence range were the Taguchi noise or environmental variables over 
which each optimal array had to be robust. A detailed description of the robust methodology used in the present 

study appears in Anderson and Keller/ 4 ^ and is termed the “Lower Order” method, while a lengthy comparison 

between the “Lower Order”, Taguchi and an alternative “High Order” method appears in Anderson and Keller/ 5 ^ 

Maximum Performance Mission - To obtain the “Optimal Adaptive” Maximum Performance 
optimal array designs, the inlet duct losses: 

Y itj = (1 -PFAVE). . (3) 

where minimized at each of the N M values of inlet throat Mach numbers and each of the N a angles-of-incidence 

to obtained the optimal array corresponding to that inlet operating condition. In a similar manner, the “Optimal 
Robust” Maximum Performance array design was determined through a search process to locate that array geom- 
etry that minimized the decision parameter: 
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N m N a 

Y M, a=FL^Ed-™)u 

M i = 1 a j=l 

where N M is the number of throat Mach number conditions in the set from Mt = 0.30 to 0.70, and N a is the 

number angles-of-incidence in the set from a = 0.0 to 20.0°. Both the “Optimal Adaptive” and “Optimal Robust” 
were both subject to the engine face distortion constraint that 

DC60 < 0.10 (5) 

and the individual Fourier harmonic 1 /2-amplitudes of distortion constraint to 

Fk 

y < 0.015 (6) 

The Optimal Robust Maximum array design was determined through a search process to have factor values of n = 
22, h = 1.95 mm, and c = 72.0 mm. In order to validate the DOE prediction results, a set of nine cases was run 
using the “Optimal Robust” array design determined from the search procedure described. The Maximum Perfor- 
mance validation cases® were organized as a full factorial array with the two mission variables at three levels 
each, i.e. 3 2 experiments. 

Comparisons between the performance results of the “Optimal Robust” and “Optimal Adaptive” 
arrays for the Maximum Performance inlet mission are shown in Figure (3). Also presented in Figures (3) is the 
inlet baseline performance at an inlet throat Mach number of 0.70 for each response. It is apparent from Figure 
(3) that flow control was able to increase total pressure recovery substantially above the baseline flow at 0.70 inlet 

throat Mach number. This was not the case for the high strength micro-vane effector units.®"® The “Optimal 
Robust” and “Optimal Adaptive” arrays provided essentially the same performance over the inlet throat Mach 

number range for 0.30 to 0.70 and angle-of-incidence range of 0° to 20°. This is not surprising, since there exists 
experimental data®) that demonstrate that a fixed secondary flow control array optimally designed can provide 
essentially the same low DC60 distortion level, i.e. DC 60 < 0.10 , over a substantial angle-of-incidence range. 
Secondary flow control in inlets is inherently robust, provided it is optimally designed. In addition, “Optimal 
Robust” and “Optimal Adaptive” arrays reduced all the Fourier harmonic 1/2-amplitudes to a value of 0.01 or 
below, which is extremely low. Although a correlation between engine face distortion and the Fourier harmonic 1/ 
2-amplitudes can not be established, reducing the Fourier harmonic 1/2-amplitudes is essentially the same at 
reducing the engine face circumferential distortion. 

Maximum Engine Stability Mission - To obtain the “Optimal Adaptive” Maximum Engine 
Stability array designs, a search was made over the factor variable space to locate that array geometry that mini- 
mized the decision parameter: 

Y.j = (DC60) i j (7 ) 

at each throat of the N M inlet throat Mach numbers and each of the N a inlet angles-of-incidence. :In a similar 

manner, the “Optimal Robust” Maximum Engine Stability array design was determined through a search process 
to locate that array geometry that minimized the decision parameter: 



where N M is the number of throat Mach number conditions in the set from Mt = 0.30 to 0.70, and N a is the 

number of angles-of-incidence in the set from a = 0.0 to 20.0°. Both the “Optimal Adaptive” and “Optimal 
Robust” arrays were both subject to the constraint that 

Fk 

y< 0.015 (9) 

while no constraint was placed on the inlet total pressure recovery (PFAVE).The Optimal Robust Engine Stabil- 
ity array design was determined through a search process to have factor values of n = 24, h = 2.0 mm, and c = 
70.2 mm. In order to validate the DOE prediction results, a set of nine cases was run using the “Optimal Robust” 

array design determined from the search procedure described. These validation cases® were organized as a full 

factorial array with the two mission variables at three levels each, i.e. 3 experiments. 

Maximum HCF Life Expectancy Mission -The “Optimal Adaptive” Maximum HCF Life 
Expectancy MSFC array was determined through a search process over the factor variable space to locate that 
array geometry that minimized the mean of the first five Fourier harmonic 1/2-amplitudes of engine face distor- 
tion, i.e 
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at each of the N M inlet throat Mach numbers and each of the N a inlet angles-of-incidence. The “Optimal 

Robust” Maximum HCF Life Expectancy MSFC array was also determined through a search process over the 
factor variable space to locate that array geometry that minimized the decision parameter: 


AT 


N„ 


M, a 


1 V 1 V 1 V ( Fk 

~~ N M " N " 5 " v 2 

M i= 1 a j = 1 k = 1 


(ID 


hj 


where N M is the number of throat Mach number conditions in the set from Mt = 0.30 to 0.70 
number angles-of-incidence in the set from a = 0.0 to 20.0°. Both the “Optimal Adaptive” and 
arrays designs were both subject to the engine face distortion constraint that 

DC60 < 0.10 (12) 

and each Fourier harmonic 1/2-amplitude satisfy the expression: 

Fk 

y< 0.015 (13) 


and N a is the 
Optimal Robust” 


The Optimal Robust HCF Life Expectancy array design was determined through a search process to have factor 
values of n = 19, h = 1.70 mm, and c = 63.0 mm. In order to validate the DOE prediction results, a set of nine 
cases was run using the “Optimal Robust” array design determined from the search procedure described. The 

Maximum HCF Life Expectancy validation cases*® were likewise organized as a full factorial array with the two 
mission variables at three levels each, i.e. 3 experiments. 


Statistical Comparison of CFD Analysis and DOE Predictions 


Presented in Anderson, Baust, and Agrell.® are the results of the statistical comparison 
between the CDF analysis and DOE predictions for the Maximum Performance, Maximum Engine Stability, and 
Maximum HCF Life Expectancy mission validation cases.These comparisons were made over the range of throat 

Mach Numbers from 0.30 to 0.70 and inlet angle-of-incidences from 0° to 20.0.° In general, the number of inci- 
dences when the comparisons were statistically different was somewhat above 5%, which is remarkably good. 
All the cases in which a statistical difference were indicated involved in the evaluation of the Fourier harmonic 1/ 
2-amplitudes of distortion. In these particular cases, the differences between the CFD analysis and DOE predic- 
tion were to too small to be of practical significance, i.e. they could not be experimentally measured. This indi- 
cates that the DOE prediction results were not statistically different from the CFD analysis results (i.e. the CFD 
analysis predictions fell within the 95% confidence interval of the DOE performance predictions). It also indi- 
cated that the optimal arrays determined by the DOE models were a statistically valid optima when compared to 
the actual CFD array analyses. The accuracy of the response surfaces determined from the DOE analysis was 
therefore more than adequate for use in determining an array optimum and for conceptual studies on the inlet- 
engine flow management system. 

Comparison of the Optimal Robust Array Designs 

Comparison of the performance of the three “Optimal Robust” array designs, i.e. the Maximum 
Performance, Maximum Engine Stability, and Maximum HCF Life Expectancy mission designs, are shown in 

Figures (4) through (8) at a throat Mach numbers of 0.30, 0.50, and 0.70 and an inlet angles-of-incidence of 0.0,° 

10.0° and 20.0°. These figures also include the baseline inlet performance, i.e. the performance of the redsigned 
M2 129 inlet S-duct without flow control. The low strength effector units used in these designs achieved a sub- 
stantial improvement in inlet total pressure recovery (PFAVE) over the baseline performance. See Figure (9). This 
differs from the performance of the high strength effector units which never increased the inlet total pressure 

recovery above the baseline value*®. Excellent engine face distortion characteristics were also achieved with the 
low strength effector units as shown in Figure (10). Although very low engine face distortion was also achieved 

with the high strength effector units®, the overall array reductions were substantially greater for the present 
designs. Presented in Figure (11) is a comparison of the first five Fourier harmonic 1/2-amplitudes for the three 
“Optimal Robust” array designs with the baseline inlet characteristics. Minimizing the mean of the first five Fou- 
rier harmonic 1/2-amplitudes resulted in a substantial reduction in the amplitudes of the first three harmonics 1/2- 
amplitudes, and very low amplitudes for the fourth and fifth harmonic components. 

By visually comparing the performance of the three “Optimal Robust” arrays designs presented 
in Figures (10) through (12), it is obvious that they are remarkably similar. This similarity can be established 
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objectively by a statistical comparison between the optimal CFD performance validations presented in Anderson, 

Baust, and Agrell.® Since each of the three (3) sets of nine (9) CFD validation cases listed in these tables were 
run at the same conditions of throat Mach number and inlet angle-of-incidence, they represent three (3) sets of 
nine (9) matched pairs of CFD observations. The data in these tables have been ordered such that the differences 
in each of the nine (9) matched pairs of CFD observations can be tested as a paired t-test. In a paired t-test, the 
mean of the sample difference and the standard deviation of the sample difference is calculated and the following 
t-statistic determined: 


where 


t* = 


1 

N 


I*,- 


L=L 



(AyAf 
(N- 1 ) 


(14) 


(Fj - F 2 ) • I s the difference of each of the N-pairs of the response variable in the two data sets. 
In comparing the mean values from two data sets (1) and (2), if the expression 


t* < t( 0.975, v p ) (15) 

is valid, the response values from the first data set are not statistically different from the response values from the 
second data set at the 95% confidence level. Likewise, the response values from the first data set are statistically 
different from the response values from the second data set at the 95% confidence level if the expression 


t* > t( 0.975, V p ) (16) 


holds. The term v is the “pooled” degrees of freedom which, for a paired t-test, is given by the expression: 

Vp = (N- 1) (17) 

and t( 0.975, V p ) is the is the 95% confidence t-value for V p degrees of freedom. 

The results of the paired t-test for the CFD validation cases are presented in Tables (10), (1 1) and 
(12) These evaluations have been organized as three sets of comparisons for two “Optimal Robust” mission array 
designs. In the first comparison, the Maximum Performance (data set 1) and Maximum Engine Stability (data set 
2) mission cases are evaluated for the seven response variables. These comparison are presented in Tables (10). In 
the second comparison the Maximum Performance (data set 1) and Maximum HCF Life Expectancy (data set 3) 
mission cases are evaluated, again for the same seven response variables. These comparisons are presented in 
Tables (1 1). Table (12) present the results for the third set of comparisons, i.e. between the Maximum Engine Sta- 
bility (data set 2) and Maximum HCF Life Expectancy (data set 3) mission cases for the same seven response 
variables. The seven response variables evaluated were the inlet total pressure recover (PFAVE), the engine face 
distortion (DC60), the first five Fourier harmonic 1/2-amplitudes of distortion (Fl/2, F2/2, F3/2, F4/2, F5/2). The 
results of this study indicated that there were no statistical significant differences between the three sets of CFD 
validations cases at the 95% confidence level. Even though there are differences in the factor variables that define 
the “Optimal Robust” array designs, these factor differences did not translate into statistically significant inlet 
performance differences over the range of throat Mach Numbers from 0.30 to 0.70 and inlet angle-of-incidences 

from 0° to 20.0°. 


CONCLUSIONS 

The fundamental importance of Genichi Taguchi’s contribution to RSM over traditional design 
approaches lies in the idea that process and product sensitivity to their environment can be incorporated into the 
optimal statistical Design-of-Experiment and subsequent analysis of the data. The Taguchi noise factors that 
cause variability in industrial design are often the environmental variables, such as temperature and humidity, 
properties of the material, and product aging. In aerodynamic design, the Taguchi noise factors can be identified 
with the mission variables, since they produce variation from the design condition. Being able to include the mis- 
sion variables directly into the inlet design process represents a major breakthrough in the area of aerodynamic 
design of inlets. The inlet system can now be designed to operate with optimal performance over a range of spec- 
ified mission variables.Taguchi’s Robust Parameter Design method, however, may not be optimal in the design 
of secondary flow arrays for inlet systems because: (a) it loses information vital to the aerodynamicist and, (b) it 
is costly. Fortunately, the important aspects surrounding Taguchi’s approach to Robust Parameter Design can 
and have been incorporated into an alternate economical approach and adapted to the inlet design problem. This 
alternate inlet design method, using a combined array approach to economical Robust Optimization Design, had a 
significant run size savings over a traditional Taguchi approach, i.e. 27 CFD experiments as compared to 45 CFD 
experiments. The combined array DOE format, in which the factor (design) variables are included with the envi- 
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ronmental (mission) variables, allows for conceptual studies to me made on the inlet-engine control system to 
determine the most efficient and cost effective management system prior to any experimentation. These concep- 
tual studies on the inlet-engine flow management system can not be made using Taguchi’s Robust Parameter 
Design methodology. 

To illustrate the potential of Robust Optimization methodology, three different mission strategies 
were considered for the subject inlet, namely (1) Maximum Performance, (2) Maximum Engine Stability, and (3) 
Maximum HCF Life Expectancy. The Maximum Performance mission minimized the inlet total pressure losses, 
the Maximum Engine Stability mission minimized the engine face distortion (DC60), while the Maximum HCF 
Life Expectancy mission minimized the mean of the first five Fourier harmonic amplitudes, i.e. “collectively” 
reduced all the harmonic 1/2-amplitudes of engine face distortion. Although the three “Optimal Robust” array 
designs were generated from three very different mission strategies, the performance achieved by these array 
designs were not statistically significantly different over the entire mission variable range. Hence one can draw 
overall conclusions with regard to microscale secondary flow control array design. Since each of the “Optimal 
Robust” array designs achieved near uniform engine face circumferential distortion, one can conclude that this 
condition represents the most robust operating state of the inlet In addition, although there is no established corre- 
lation between the seven response variables tested in this study, it is clear that numerically optimizing the second- 
ary flow array designs over the mission variable range drives all the response variables to their optimal levels. In 
Robust Optimization, therefore, the system is greatly simplified because the inlet response variables are collec- 
tively managed. Had the array not been numerically optimized, the response variables would have to have been 
individually managed, and the inlet management system would become far more complex. 

In general, the performance differences between the “Optimal Adaptive” and “Optimal Robust” 
array designs were found to be marginal. This suggests, that “Optimal Robust” open loop array designs can be 
very competitive with “Optimal Adaptive” close loop designs. Thus the design concept of Robust Optimization, 
where the variability cause by the mission parameters is minimized, is achievable and valid. This design concepts 
results in a simpler inlet management systems as compared to closed-loop active feedback control. Under closed- 
loop active feedback control, the state of the inlet must be sense and a corresponding reaction to that state must be 
implemented for each response variable. Taking advantage of the Optimal Robustness properties of microscale 
secondary flow control therefore can lead to simpler, more reliable and more cost effect inlet management sys- 
tems. 
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Figure (1): Location of effector region within inlet 
S-duct configuration 


Factor 

Range 

Number of Vane Effectors, n 

13 to 27 

Effector Vane Height (mm), h 

1.0 to 2.0 

Effector Chord Length (mm), c 

36.0 to 72.0 

Inlet Throat Mach Number, Mt 

0.30 to 0.70 

Inlet Angle-of-Incidence (degs.), oc 

0.0 to 20.0 


Table (1): Factor variables which establish the DOE 
design matrix. 


Variable 

Value 

Effector Vane Thickness (mm), t 

0.138 

Vane Angle-of-Incidence (degs), (3 

5.0 

Inlet Total Pressure (lbs/ft 2 ), Pt 

10506.0 

Inlet Total Temperature (°R), Tt 

517.0 

Inlet Angle-of-Yaw (degs), y 

0.0 


Table (2): Variables held constant. 


Response 

Nomenclature 

Engine Face Total Pressure Recovery 

PFAVE 

Engine Face Distortion 

DC60 

1st Fourier Harmonic 1/2-Amplitude 

Fl/2 

2nd Fourier Harmonic 1/2- Amplitude 

F2/2 

3rd Fourier Harmonic 1/2-Amplitude 

F3/2 

4th Fourier Harmonic 1/2- Amplitude 

F4/2 

5th Fourier Harmonic 1/2- Amplitude 

F5/2 


Table (3): DOE response variables 


Config. 

n 

h 

c 

Mt 

a 

nvg701 

13 

1.0 

36.0 

0.30 

20.0 

nvg702 

27 

1.0 

36.0 

0.30 

0.0 

nvg703 

13 

2.0 

36.0 

0.30 

0.0 

nvg704 

27 

2.0 

36.0 

0.30 

20.0 

nvg705 

13 

1.0 

72.0 

0.30 

0.0 

nvg706 

27 

1.0 

72.0 

0.30 

20.0 

nvg707 

13 

2.0 

72.0 

0.30 

20.0 

nvg708 

27 

2.0 

72.0 

0.30 

0.0 

nvg709 

13 

1.0 

36.0 

0.70 

0.0 

nvg710 

27 

1.0 

36.0 

0.70 

20.0 

nvg711 

13 

2.0 

36.0 

0.70 

20.0 

nvg712 

27 

2.0 

36.0 

0.70 

0.0 

nvg713 

13 

1.0 

72.0 

0.70 

20.0 

nvg714 

27 

1.0 

72.0 

0.70 

0.0 

nvg715 

13 

2.0 

72.0 

0.70 

0.0 

nvg716 

27 

2.0 

72.0 

0.70 

20.0 

nvg717 

13 

1.5 

54.0 

0.50 

10.0 

nvg718 

27 

■a 

54.0 

0.50 

10.0 

nvg719 

20 

1.0 

54.0 

0.50 

10.0 

nvg720 

20 

2.0 

54.0 

0.50 

10.0 

nvg721 

20 

1.5 

36.0 

0.50 

10.0 

nvg722 

20 

1.5 

72.0 

0.50 

10.0 

nvg723 

20 

1.5 

54.0 

0.30 

10.0 

nvg724 

20 

1.5 

54.0 

0.70 

10.0 

nvg725 

20 

■a 

54.0 

0.50 

0.0 

nvg726 

20 

1.5 

54.0 

0.50 

20.0 

nvg727 

20 

1.5 

54.0 

0.50 

10.0 


Table (4): “Lower Order” Central Composite 
Face-Centered (CCF) combined array DOE format 
involving factor (design) variables and environ- 
mental (mission) variables 
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a=0.0° a=10.0° a=20.0° 









Figure (2): Graphical representation of the “Lower Order” Central Composite Face-Cen- 
tered (CCF) combined array DOE format involving factor (design) variables and environ- 
mental (mission) variables. 
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Baseline Perf., Mt = 0.70 

Optimal Robust, Mt = 0.30 

Optimal Robust, Mt = 0.50 

Optimal Robust, Mt = 0.70 

° Optimal Adaptive, Mt = 0.30 

^ Optimal Adaptive, Mt = 0.50 

^ Optimal Adaptive, Mt = 0.70 



Inlet Angle-of- Incidence (degs.), cl 


(a) Engine Face Total Pressure Recovery 



Inlet Angle-of-Incidence (degs.), a 

(b) Engine Face DC60 Distortion 



(c) 1st Fourier Harmonic 1/2 Amplitude 



Inlet Angle-of-Incidence (degs.), a 

(d) 2nd Fourier Harmonic 1/2 Amplitude 
Figure (3): Comparison of “Optimal Robust” and 
“Optimal Adaptive” Maximum Performance inlet 
mission array designs 


Baseline Perf., Mt = 0.70 
Optimal Robust, Mt = 0.30 
Optimal Robust, Mt = 0.50 


Optimal Robust, Mt = 0.70 

° Optimal Adaptive, Mt = 0.30 
^ Optimal Adaptive, Mt = 0.50 



Inlet Angle-of-Incidence (degs.), oc 


(e) 3rdFourier Harmonic 1/2 Amplitude 



Inlet Angle-of-Incidence (degs.), a 

(c) 4th Fourier Harmonic 1/2 Amplitude 



Inlet Angle-of-Incidence (degs.), a 

(g) 5th Fourier Harmonic 1/2 Amplitude 

rt 

fa 0.050 

oT 

T3 

a 0.040 

a 

£ 

< 0.030 

rj 

0.020 

0 

1 0.010 

* > 

a 0.000 

I 0 2 4 6 8 10 12 14 16 18 20 

Inlet Angle-of-Incidence (degs.), a 

(h) Mean Fourier Harmonic 1/2 Amplitude 
Figure (3): Comparison of “Optimal Robust” and 
“Optimal Adaptive” Maximum Performance inlet 
mission array designs, continued. 
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(b) Optimal Robust Max. Performance 
Figure (4): Near wall inlet S-duct streamline traces, 
Mt = 0.30, a = 0.0°. 




(b) Optimal Robust Max. HCF Life Expct. 
Figure (5): Near wall inlet S-duct streamline traces, 
Mt = 0.70, a = 20.0°. 



(a) Baseline Inlet (b) Optimal Robust 

Figure (6): Total pressure recovery contours, Opti- 
mal Robust Maximum Performance array, Mt = 

0.30, a = 0.0°. 



(a) Baseline Inlet (b) Optimal Robust 


Figure (7): Total pressure recovery contours, Opti- 
mal Robust Maximum Engine Stability array, Mt = 

0.50, a = 10.0°. 



(a) Baseline Inlet (b) Optimal Robust 


Figure (8): Total pressure recovery contours, Opti- 
mal Robust Maximum HCF Life Expectancy array, 

Mt = 0.70, a = 20.0°. 
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Response 

MEAN of 
Sample 
Diff. 

t* 

t(0.975,8) 

Comment 

PFAVE 

0.0 

0.0 

2.306 

Not Diff. 

DC60 

0.0046 

0.1734 

2.306 

Not Diff. 

Fl/2 

-0.0008 

0.3916 

2.306 

Not Diff. 

F2/2 

-0.0011 

0.4242 

2.306 

Not Diff. 

F3/2 

-0.0006 

0.6299 

2.306 

Not Diff. 

F4/2 

0.0004 

0.8433 

2.306 

Not Diff. 

F5/2 

0.0 

0.0 

2.306 

Not Diff. 


Figure (9): Effect of “Optimal Robust” array 
designs on inlet pressure recovery (PFAVE), Mt = 

0.70, a = 20.0°. 


Table (5): Paired t-test statistical evaluation 
between the (9) optimal Max. Performance cases 
and the (9) Max. Engine Stability cases 



Baseline Max.Perf. Max. Stability Max. HCF Life 


Figure (10): Effect of “Optimal Robust” array 
designs on engine face distortion (DC60), Mt = 

0.70, a = 20.0°. 


Response 

MEAN of 
Sample 
Diff. 

t* 

t(0.975,8) 

Comment 

PFAVE 

0.0008 

0.3190 

2.306 

Not Diff. 

DC60 

0.0043 

0.3911 

2.306 

Not Diff. 

Fl/2 

0.0008 

0.2600 

2.306 

Not Diff. 

F2/2 

-0.0006 

0.2611 

2.306 

Not Diff. 

F3/2 

-0.0006 

0.4167 

2.306 

Not Diff. 

F4/2 

0.0001 

0.1421 

2.306 

Not Diff. 

F5/2 

0.0002 

0.2667 

2.306 

Not Diff. 


Table (6): Paired t-test statistical evaluation 
between the (9) Max. Performance cases and the (9) 
Max. HCF Life Expectancy cases. 


0.05 

0, 0.04 

| 

a 

S 0.03 
CJ 

I 0.02 
© 

X 0.01 
0 

Fl/2 F2/2 F3/2 F4/2 F5/2 FM/2 










r 

.*lJ 

flwHiil m 


□ Baseline 

■ Max. Perf. 

■ Max. Stability 

□ Max. HCF Life 


Harmonic Component 


Figure (11): Effect of “Optimal Robust” array 
designs on the Fourier harmonic 1/2 amplitude 

(Fk/2), Mt = 0.70, a = 20.0°. 


Response 

MEAN of 
Sample 
Diff. 

t* 

t(0.975,8) 

Comment 

PFAVE 

0.0008 

0.3590 

2.306 

Not Diff. 

DC60 

0.0084 

0.8685 

2.306 

Not Diff. 

Fl/2 

0.0156 

0.4794 

2.306 

Not Diff. 

F2/2 

0.0006 

0.7647 

2.306 

Not Diff. 

F3/2 

0.0 

0.0 

2.306 

Not Diff. 

F4/2 

-0.0003 

0.3333 

2.306 

Not Diff. 

F5/2 

0.0002 

0.0544 

2.306 

Not Diff. 


Table (7): Paired t-test statistical evaluation 
between the (9) Max. Engine Stability cases and the 
(9) Max. HCF Life expectancy cases. 
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